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Abstract. A comparison of the recent age and temperature estimates of 
the Geminga neutron star with cooling models is presented. This star is 
already in the photon cooling era and it is shown that its temperature can 
be understood within both the slow and fast neutrino emission scenarios 
and consequently will not allow discrimination between these two scenarios. 
However in both cases agreement of the theoretical cooling curves with the 
observed temperature depends crucially on the presence of nucleon padring 
in most, if not all, of the core. 


1. Introduction 

The study of the thermal evolution of young neutron stars offers the possibility to 
obtain unique information about the structure of compressed cold nuclear matter. The 
early cooling after the supernova explosion is driven by neutrinos whose emission rate 
is a very sensitive function of the state of that matter. Moreover, both the neutrino and 
the later photon coolings are strongly affected by the occurence of pairing (a la BCS) 
of the baryonic components in the star’s core. These two effects give us two handles to 
extract information about nuclear matter through comparison with observed neutron 
stars of known age. The announcement of the definitive identification of Geminga as 
a spinning magnetized neutron star by Halpem & Holt this spring gave rise to great 
excitement in the astrophysical community. This object, which has been a mystery for 
many years x , is suddenly revealing a wealth of informations thanks the ROSAT and 
GRO observatories. We present here a comparison of the Geminga observation with 
our latest cooling models and show what can be learned from it. These results will 
hopefully stimulate more detailed works on the various aspects of the models presented. 

1 Geminga stands for Gemini Gamma-ray source , but it also means ‘does not exits’ or ‘is not there’ 
in Milanese dialect [1], 
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2. Geminga 


2.1. The story 

Geminga was discovered in 1972 by the 7 -ray telescope on SAS-2 [2] and was observed 
again five times by COS-B between 1975 and 1982. It was the second bnghtest source 
in the sky at energies above 100 MeV in the second COS-B catalog [3] where t was 
labelled 2CG195+04. Gamma-ray detectors did not have enoug re ®° u 

to allow for search of an optical counterpart, but the situation changed when the 
EINSTEIN satellite could find four soft X-ray sources within t e r»y error ojo 
A ndvsis showed that only one of these sources, 1E0630+178, could be the counterpart 
Arr'l 95+04 111 The angular resolution of EINSTEIN’s HRI (90% confidence radars 
SSK op« search. 1E0630+178 points to a blank field on the 
Palomar Sky Survey plate but deeper observations revealed two stars m the error 
circle and the faintest one, with r-magnitude ~ 25.5, is presently considered as the 
optical counterpart of Geminga [4,5). It thus took ten years to find G< ^ “ S °* 

X rays and more than fifteen years to identify it in the optical range. It has not yet 

been detected at radio wavelength. 

Early on it was suspected that Geminga was an isolated neutron star but other 
models had been proposed. The situation dramatically changed this spnng when 
Halpera k Holt [6] discovered pulsations at a period of 0.237 sec in their soft X 
ray ROSAT observation of Geminga. Using this result, Bertsch et si [7 found the 
same period in the 1991 GRO 7-ray observation 

it in the archival COS-B data. This definitely identified 2CG195+04 and 1K) 630 + 
as the same object and proved unequivocally that Geminga is an isolated spinning 

magnetized neutron star. 

2.2. The data 

2 2 1 Geminga ’* age. The Geminga period as obtained in [6] is P = 0.2370974 s ± 
01 Is. The earlier observations by COS-B and later by GRO have shghtly different 
P’s which lie very accurately on straight line (see [8]) and give a practically constant 
period derivative P = 1.0977 x 10- s s- |9] over a span of 16 yearn. From these we 
obtain its ‘spin-down age’ r = &P/P = 3.4 x 10 s years with n = 3 
for magnetic dipole breaking). The relation between the spin-down age and the real 
age of f pulsar is a delicate problem, refer to [10,11,12] for discussions. We will take 
as extreme values the two case of n = 2 and n = 4, giving 

2.3 x 10 5 yrs < t < 6.8 x 10 5 yrs, C 1 ) 

-the upper value being probably an overestimate. 

2.2.2. Geminga’s temperature. In the first estimate of the star temperature Halpem 
& Holt [6] fitted the two apparent components of the spectrum by a black y an 
a power law spectra. The best fit gave a blackbody temperature of T - 3 - 4 X 1U 
K In a second, more detailed, analysis Halpem and Ruderman (13] replace the power 
law component by a second black body with higher temperature and argue that this 
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hotter component (T = 3x 10 6 K) is due to emission from a reheated polar cap (this 
temperature is too high to be explained only by anisotropic magnetic heat transport). 
They obtain for the main surface emission a temperature of T = (5.2 ± 1.0) X 10 K. 
From the ratio of the luminosities of these two components they conclude that the 
ratio of the areas of the hot and cold emitting regions is about 3 x 10 5 . Moreover it 
is likely that the surface is not emitting uniformly and that some colder region is not 
been seen, making the concept of ‘surface temperature’ an ambiguous one. However 
cooling calculations give as an output the effective temperature, i.e, a total luminosity, 
and the presence of a cooler region would reduce the luminosity and make the effective 
temperature somewhat lower than the above value of (5.2 ± 1.0) x 10 s K. 


All these analyses use black-body spectra but the surface of a neutron star cannot 
be expected to be a perfect black body. Romani [14], Miller [15] and Shibanov et al 
[16] have calculated more realistic spectra for various surface chemical compositions 
with ([15,16]) or without ([14]) magnetic field. The general trend of these results, for 
H, He or Fe atmospheres, is that there is an excess emission in the Wien tail of the 
spectrum compared to the blackbody which falls within the EINSTEIN and ROSAT 
detector ranges. This excess is reduced if metals are present (due to absorption edges) 
or when the effect of the magnetic field is taken into account. Use of these spectra 
will lower the measured temperature. Finally some contamination from a surrounding 
nebula and/or surface reheating by gamma-rays or particles from the magnetosphere 
cannot be excluded. 

We will take for comparison with our calculations an effective temperature of 
4 - 6 x 10 3 * 5 K (= lO 5 * 6 - 5 ’ 75 K) but insist that it has to be taken as an upper limit. 


3. Neutron star cooling 

After being formed hot in a supernova explosion, a neutron star cools by neutrino 
emission from its core and crust and by photon emission from its surface. Neutrino 
emissivities are proportional to Tf or Tf, while photon emissivity is proportional to 
2? ~ T?* (2* is the interior temperature and T, the surface temperature; T # ~ 

[17]), consequently photons will eventually dominate when T{ has dropped sufficiently. 
The shift from neutrino to photon cooling usually occurs at about 10 5 years of age. 
Geminga is in the photon cooling era. 

8.1. Neutrino processes 

The dominant neutrino emission processes occur in the core of the star and are variant 
of beta and inverse beta decays. Table 1 shows the approximate relevant emissivities 
for c ompa rison - One can divide them into slow and fast neutrino emission, depen din g if 
they involve four or two baryons. The difference between fast and slow emission comes 
mainly from phase space considerations, and the precise values depend on form factors, 
strangeness violation and strong interaction corrections. We refer to the recent review 
by Pethick [18] for a detailed discussion. There is still doubt about which processes 
can act ually occur, but the number of proposed channels for fast cooling is now so 
large that it is becomming difficult to believe that none of these is permitted and that 
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Table 1. Some core neutrino emission processes and their emissivities. Hy- 
perons, if present, can go into either modified or direct Urea processes with 
emissivities slightly lower than the corresponding purely nucleonic processes 
[19]. T 9 is the temperature in units of 10 9 kelvins. 


Process Name 


Process 

Emissivity Q v 
(erg/sec/cm 3 ) 

Modified URCA 

f n + n'— mi' + p + e + v t 
1 n' + p + e~ — * n' + n + u e 

10 20 • 

T* 

[20] 

K-condensate 

: 

+ K " n + e“ + v e 

+ e" -*• n + K ~ + v e 

- 10 24 • 

It 

[21] 

7r - condensate 

1 n + i — *n + e +i/ e 
1 n + e" -* n + tt _ + v c 

^ 10 26 • 

7? 

[22] 

Direct URCA 

i 

[ n —* p + e~ + 

1 p + e~ —* n + u e 

~ 10 27 - 

TS 

[23] 

Quark URCA 

< 

) d — » u + e" + Pe 
( u + e~ d + v e 

~ 10 26 q 

■X 

[24] 


neutron star cooling follows the old ‘standard model’ with only the modified Urea 
process. Awaiting for a definite argument on this point we will still consider both the 
fast and slow cooling scenarios. 

Deconfined quarks may be present in the center of massive neutron stars and 
are also copious neutrino emitters. They thus belong to the fast cooling scenario but 
obviously require a separate treatment. We will not consider them explicitely here. 

S.2. Nucleon pairing 

Pairing of nucleons in the neutron star core has a dramatic effect on the cooling, both 
by suppressing the neutrino emission and the specific heat. It was shown in [34] that 
the temperature at ages between lO’-lO 5 years of a fast cooling neutron star is actually 
determined by the temperature T c at which its core nucleons are paired. Since T c is 
density dependent it is its lowest value (usually in the very center of the star) which 
is relevant. Theoretical calculations of T c are extremelly difficult and the presently 
published values are still very uncertain. In the core, the protons are paired in the 
1 5o partial wave while the neutrons are in the 3 P 2 partial wave. x 5o pairing is easier to 
study than 3 P 2 , but unfortunately, in the case of protons, it depends strongly on the 
relative concentrations of protons and neutrons since protons see mostly neutrons. For 
neutron 3 P 2 pairing the proton fraction x v is not so important as long as it is low, but 
s interaction is much more difficult to handle than the x So interaction (moreover 
the tensor interaction also introduces an important 3 F 2 contribution). 

Figure la shows the results of calculations of proton pairing critical temperatures. 
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a) Proton % pairing « Neutron % pairing 



Figure 1. a) Proton J 5 0 pairing critical temperatures. CCY : [25], 
T73 : [26], NS : [27], AO : [28], WAP : [29]. 

b) Neutron pairing critical temperatures. HGRR : [30], T72 . 
[31], AO : [33]. The two dashed curves show the results of T72 and AO 
when the effective mass is fixed to the free mass value. 


One sees that their maxima differ by a factor five and that there is also a large spread in 
the Fermi momentum at which T c vanishes. The most recent calculation [29] is the first 
one to take properly into account the presence of the neutrons; it gives a lower T c than 
the previous ones but which does not decrease with density as fast as in the previous 
^/•ulatinng - However this results may c han ge when a different proton fraction is used 
(T. Ainsworth, private communication). 

Figure lb shows the three calculations of neutron critical temperatures pub- 
lished to date. One sees that the maximum value of T e as well as the Fermi momentum 
at which T c vanishes differ enormously from one calculation to the other. The curves 
T72 [31] and AO [33] use the same nuclear potential but different many-body methods 
and give the sam e maximum T c but very different vanishing densities. A fourth calcu- 
lation [32] of the pairing interaction which includes medium effects not considered in 
the other ones seems to give a very high T e and a strong dependence of the vanishing 
density on the potential. From this one sees that the value of T c for neutron Pj pairing 
and the density range where pairing occurs are still open questions. 

If hyperons are present one can expect that they also pair for the same reasons 
as the nucleons do. With regard to quarks, pairing is also very probable [35]. 

The effect of pairing is to reduce the phase space available for excitations. As a 
result both the specific heat of the paired component and the neutrino processes to 
which it contributes will be suppressed. We treat this suppression by simply multi- 
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Figure 2. Cooling by the direct Urea process and neutron 
3 _P 2 pairing suppression. The dashed curve show the cooling by direct 
Urea with no core pairing. All continuous curve have direct Urea and 
core neutron pairing with various T e from Figure lb); the 0.3HGRR and 
0.1HGRR have pairing from HGRR with T e multiplied by 0.3 and 0.1 re- 
spectively. The dotted curve show the cooling with only the modified Urea 
and no core pairing for comparison. In all cases the crust neutrons are paired 
with T c from [36]. The cross shows the Geminga value, the other three data 
points are as in [34] (they should also be plotted with an uncertainty on the 
age). 


plying the corresponding normal values by a Boltz mann factor exp(— A/fcjgT). This 
correct when T 4HT e but overestimate the suppression near T e since there the gap 
vanishes. Accurate suppression factors are bring calculated by Levenfish & Yakovlev 
and will be available soon for a more quantitative treatment. - 


4. Comparison with Geminga 
4.1. Fast cooling scenario 

In the fast cooling models one accepts early enhanced neutrino emission, by the direct 
URCA process or by some ‘exotic’ matter, but must then suppress the emi ssion by 
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Table 2. Some properties of the EOS’s used. Columns two to four list the 
mass, central density and central proton fraction of a maximum mass star. 
Columns five to seven give properties of a 1.4M 0 star : radius, central density 
and central proton fraction. (The PAL EOS’s are labelled PALij, i,j— 1,2,3, 
where i refers to the symmetry energy function and j to the compression 
modulus) 


EOS 

Maximum mass star 
•Mmax(A/0) Pmax(fm ) >(%) 

R(km) 

1.4 Mq star 

Pc(fm~ 3 ) 


FP [38] 

1.79 

1.18 

™ 0 

10.85 

0.69 

2.0 

WFF(avl4) [39] 

2.10 

1.25 

4.8 

10.60 

0.64 

9.6 

MPA [40] 

2.44 

0.89 

19.0 

12.45 

0.41 

9.0 

PAL32 [41] 

1.68 

1.51 

15.2 

11.02 

0.74 

11.4 

PAL33 [41] 

1.90 

1.24 

14.0 

11.91 

0.54 

9.7 


nucleon superfluidity (Page & Baron [37] and Page & Applegate [34]). Without su- 
perfluidity suppression the resulting surface temperature is much below the estimated 
value for any neutron star observed so far. Figure 2 show typical cooling curves, from 
[34], with the direct Urea process and various Z P^ neutron gaps. A large gap as calcu- 
lated by Hoffberg et al [30] almost completely turn off the core emission while the small 
gap from Takatsuka [31] has almost no effect because most of the core remains normal. 
Intermediate values give reasonable agreement with the Geminga temperature. See 
[34] for more details. 

The fast neutrino emission is so efficient that any region in the star which re m ains 
normal will drive the surface temperature well below the observed temperature of 
Geminga. Consequently within the fast cooling scenario one of the neutrino emitting 
components must be paired at the highest densities reached in the core. If hyperons 
are present the same conclusion holds since even a small amount of A will allow a 
direct Urea reaction [19]. If both A and £“ axe present they participate into a purely 
hyperonic direct URCA process and thus one of them must be paired. One may expect 
hyperons to have a lower T e if they have weaker interactions than nucleons and if this 
is the case what controls the cooling is the A or E“ pairing critical temperature. 

4-2. Slow cooling scenario 

While with the fast cooling scenario one needs pairing to reduce the cooling, 
within the slow cooling scenario one encounters the opposite problem : neutrino emis- 
sion is not sufficient and one needs to speed up the cooling during the photon era by 
reducing the specific heat through nucleon pairing. This slow cooling scenario, often 
called the ‘standard model’, assumes that the only core neutrino emission is by the 
modified Urea process and the two (less efficient) similar processes with neutral cur- 
rents. Comparison with the standard calculations of Nomoto and Tsuruta [42] (or 
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Table 3. Normal specific heat, at T = 10 9 K, of neutrons, protons and 
electrons in the core and crust of a 1.4M© neutron star built with the five 
EOS’s used. (Units are 10 38 ergs K -1 , i.e. C V (T ) = (Table — Entry) x 
(T/10*K) x 10 38 ergs K~ l ). ' 


EOS 

Core components 

Crust components 

• 

n 

P 

e 

n 

e 

FP 

9.00 

2.89 

0.50 

0.97 

0.025 

WFF(avl4) 

9.31 

2.81 

0.51 

0.72 

0.018 

MPA 

11.9 

3.90 

0.67 

1.60 

0.044 

PAL32 

10.5 

3.58 

0.67 

1.60 

0.044 

PAL33 

9.66 

3.41 

0.68 

1.16 

0.031 


similar resluts by Van Riper [43]) shows that some EOS’s can produce a star with the 
observed temperature of Geminga, but other still give too high a temperature. How- 
ever these authors do not consider the enormous uncertainty in the density at which 
the gap vanish : their ‘successful’ models are simply the ones with stiff EOS (i.e. low 
central density) which have hence pairing in the full core. To avoid this misleading 
interplay between the EOS stiffness and the density dependence of the gap we take 
density independent gaps, therefore forcing pairing in the full core. We consider five 
different EOS’s from modem calculations, relativistic and non-relativistic, which en- 
compass a large range of stiffness and proton fraction. Their properties are s ummari zed 
in Table 2. We reject EOS’s with proton fractions large enough to allow the direct Urea 
process in a 1.4M© star, but do consider MPA, PAL32 and PAL33 which allow it only 
at higher masses. Table 3 shows the contribution to the specific heat of the various 
components in a 1.4 M© star at a temperature T = 10 9 K without pairing. Pairing 
will suppress C v exponentially when T<T C and the corresponding specific heat will 
practically disappear. In all our calculations the crust neutrons are paired using the 
gap from [36] and their contribution to C v is thus strongly reduced; the crust electrons 
have a negligible contribution as well as the crust lattice. One can see from the table 
that the core neutrons contribute about | of the total specific heat, the protons j and 
the core electrons about 0.5 %. These fractions are suprisingly almost independent of 
the EOS. 

The cooling curves for 1.4 Af© stars built with our five EOS’s are compared in 
Figure 3 with the Geminga observation. When no core pairing occurs (a) the theoretical 
results are consistent only with the higher surface temperature T, and the older age: 
taking into account that this T, is certainly an overestimate and considering also that 
several reheating mechanisms [44,45] should raise the temperature of the theoretical 
curves at this age, one can state that Geminga’s temperature is incompatible with 
these cooling models (unless Geminga’s age is underestimated). With padring of the 
protons (b) the discrepancy increases because of the suppression of the early neutrino 
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L °3 Js (X) .. .. LogT s (K) 



Log t (years) 


Log t (years) 


Figure 3. Cooling by the modified Urea process without/ with 

pairing(s) and five different EOS’s : FP , WFF(avl4) , MPA 

.... -PAL32 — • — and PAL33 — • • — . In all cases 'the crust neutrons 
are paired with T e from [29]. The cross shows the Geminga value. 

a) No core pairing at alL 

b) Protons paired with density independent T c = 2 x 10 9 K. 

c) Core neutrons paired with density independent T e = 2 x 10 9 K. 

d) Protons and core neutrons paired with density independent T c = 2 x 10 9 
K. 
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cooling and the only small (~ 25%) decrease of the specific heat in the photon cooling 
era. When neutrons are paired (c) the reduction of C v is large enough to accomodate 
the observed temperature. If both neutrons and protons are paired (d) within the 
whole core C v is cut down by a factor 20 (only the electron contribution is left) and 
the temperature drop in the photon cooling era is extremely fast; however reheating 
mechanisms could be very efficient here, due to the very low C„, and rise the curve 
significantly. The curves for the various EOS’s axe very similar when identical T c are 
used, showing that the EOS-dependent results obtained in other works (e.g. [42,43]) 
are mostly due to the density dependence of T e that the authors choose. The only thing 
that one can say about the value of T c is that it has to be high enough for the specific 
heat to be sufficiently suppressed, the value of 2 x 10 9 K used is only illustrative. 

If one relieves the approximation of a density independent gap two things can 
happen. Close to the crust-core boundary where neutron pairing shifts from 1 So to 
z Pi there may be a thin layer where T c is vanishingly small : this layer would contributes 
to C v when the deeper parts are paired, but it comprises only a small shell. In the 
center of the star the density gradient is very small and a large mass is at density close 
to the maximum density : the neutron Fermi momentum is above 90 % of its center 
value in a region of mass larger than 0.3 M 0 for a 1.4 M e star from PAL33 and a 
mass larger than 0.6 M e with WFF(avl4). Therefore if the central density is above 
the vanishing density for neutron pairing it cannot be much above it if we want a large 
enough reduction of C v . 

The envelope calculation used does not take into account the effect of the magnetic 
field on the heat transport. It has been argued that this effect is small [46] and hence 
our results should not change significantly when the magnetic field effects are included. 
A detailed analysis is in progress and will be published later. 


5. Conclusion 

We have compared the recent temperature measurement of the Geminga neutron star 
with cooling models and found that, since this star is old enough to be in the photon 
cooling era, both fast and slow neutrino emission mechanisms can explain it. One there- 
fore cannot draw any conclusion about neutrino emission from dense nuc l e a r matter 
using this observation only. However, a crucial feature in both types of models is that 
they need nucleon padring in most, if not all, of the core. With fast neutrino cooling 
nucleon padring is needed to stop the early cooling which, without this, would produce 
a star of temperature much lower that the observed one. If the fast neutrino emission 
is from hyperonic processes it is possible that the suppression we observe is due to 
hyperon superfluidity. Since fast neutrino emission occurs down to the very center of 
the core these scenarios need pairing up to the highest density reached in this object. 
With the slow cooling model (‘standard model’) superfluidity is also needed, but for a 
different reason. The observed temperature is below what the simple model predicts, 
but since this stair is cooling by photon emission we can accelerate the cooling at this 
time by decreasing the specific heat through pairing. If we accept am age of 3 x 10 5 
years and a temperature of 5 x 10 s K then the specific heat must have been reduced to 
about 25% of its normal value. This cam be obtained either by pairing of the neutrons 
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in the whole core or by a combination of both neutron and proton pairing, but even 
in this case most of the neutrons must be paired. If both neutrons and protons are 
paired in the whole core then a substantial amount of reheating is needed, but several 
possible mechanisms have been proposed and they can provide sufficient reheating. 

As a final point it should be mentioned that this results will have important 
consequences in modeling the cooling of other neutron stars. For example the cooling 
calculations of [34] indicate that the Vela pulsar and PSR0656+14 have temperature 
slightly lower than what the standard model predicts without superfluidity: including 
some pairing in agreement with the present analysis will raise the theoretical prediction 
and make the discrepancy larger, reinforcing the case for the occurence of fast cooling 
in these two objects. 
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ABSTRACT 

Recent work has shown that many modern nuclear equations of state give P ro “ n /^.°“^^ Sj 
fnr thp direct Urea orocess to operate in the interiors of neutron stars. We include the direct Urea 
in calculations of neltron star coding and find that the surface temperature of a youngneutron star 
JJ mtastroDhicallv after ~ 10 2 yr if the direct Urea process is allowed and nucleons do not become super- 
JffHSS SKSity occurs throughout the direct Urea region, the surface temperature drops to a 
Si transition temper .tier -10' yr and decreet slowly f or the nn« 

~ 10* vr at which time surface photon cooling takes over. The surface temperatures of all of the candidates 
for thermal X-ray emission from cooling neutron stars are determined by the transition temperature for super- 
fuidh^nd alZt nothmg dse if these stars cool by direct Urea. By comparison with observational data, we 

temporaturcs of the order of 10’ K are required in the whole dtreot Ujt> arc 
core. If neutron stars are created with a range of masses, it is possible that the 

Urea and the less massive do not, resulting in two populations of stars differing in their thermal properties. 
Subject headings: dense matter — stars: neutron — X-rays: stars 


1. INTRODUCTION 

Recent work (Lattimcr et al. 1991) has shown that many 
modern nuclear equations of state have large enough sym- 
metry energies and predict sufficiently high central densities to 
allow the direct Urea process to operate in the interiors of 
neutron stars. The neutrino emissivity of the dirt* JJrca 
process is orders of magnitude larger than that of the modified 
Urea process and those in all scenarios involving nonnudeomc 
decrees of freedom (pion condensate, kaon condensate, quark 
matter etc.). The Lattimer et al. (1991) result clearly necessi- 
tates a’ complete reexamination of the theory of neutron star 
cooling. In this Letter we report the first results from our 
program to perform such a reanalysis. . 

In the absence of nucleon superfluidity, the interior of a star 
which cools by the direct Urea process cools very rapidly. The 
surface temperature of such a star drops catastrophically once 
the crust completes its thermal relaxation, which typically 
reauires ~ 10 2 yr. Neutron or proton superfluidity suppresses 
neutrino emission and halts the rapid cooling because 
nucleons must be excited above the gap in order to participate 
in the Urea process. Nucleon superfluidity stabilizes the 
surface temperature at a value determined by the superfluid 
transition temperature in the interior and almost nothing else. 
The aualitative effect of superfluidity shutting off rapid neu- 
trino emission was noted by Page & Baron (1990), who studied 
the effect of kaon condensation on neutron star cooling. 

The thermal histories of stars which cool by direct Urea and 
those which do not are very different. If neutron stars are 
created with a range of masses comparable to the mass range 
of the iron cores of their progenitors, it is possible that the 
more massive stars cool by direct Urea and the less massive 
cool by less efficient means. We illustrate the dramatic effect of 
direct Urea cooling by choosing an equation of state which has 
a minimum mass of M cr * - 1.35 M© for the onset of direct 
Urea cooling and presenting cooling curves for store in the 
mass range 1.0 M© < M < 1.7 M 0 . The difference is striking, 
but the quantitative degree to which neutron store fall into two 


thermal populations requires the careful treatment of the role 
of hyperons and a possible hyperon Urea process (Prakash et 
al. 1992). 

2. THE MODEL 

We solve the relativistic equations for heat transport 
through an unmagnetized neutron star using the Hcyney-type 
stellar evolution code described in detail in Page (1989) and 
Page & Baron (1990, 1991). We use a nuclear equation of state 
with a symmetry energy allowing the direct Urea process, and 
we use the Lattimer et al. (1991) value for the neutrino emiss- 
ivity due to the direct Urea process. Other than these, the only 
change is that we use improved calculations of the Coulomb 
logarithm for electron conductivity in the crust liquid phase by 

Yakovlev (1987). . , 

The nuclear symmetry energy is the most important factor 
determining whether or not the direct Urea process is allowed 
in the star. We parameterize the core equation of state using 
the functional form suggested by Prakash, Ainsworth, & Latti- 
mer (1988). These authors define the nuclear symmetry 
energy S(n) by writing the energy per nucelon as 

£(n, x) = £(n, 0.5) + S(«X1 - 2x) 2 , 

where x = Z/A is the proton fraction, n is the number density 
of nucleons, and £(n, 0.5) is the energy per nudeon of sym- 
metric nuclear matter. We choose a symmetry energy of S(n) - 
3<Xn/n 0 ) 07 MeV, where n 0 = 0.16 fm' 3 is the saturation 
density of symmetric nuclear matter. This choice is typical of 
the symmetry energies found in the relativistic equations of 
state analyzed by Lattimer et aL (1991). We take the other 
parameters from the second line of Table 1 in Prakash et al. 
(1988). These parameters give a compression modulus of K 0 = 
180 MeV, which is the number advocated by conventional 
wisdom, a maximum neutron star mass of 1.7 M©, and a criti- 
cal mass above which the direct Urea process is allowed of 

Man = 1-35 Af©. , L . . . 

Superfluidity alters both the specific heat of the supernuia 
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species and the neutrino emissivity because of the gap that it 
introduces into the excitation spectrum. At low density the 
pairing is predicted to be in the ’S 0 state for both neutrons and 
protons. At high density the neutron pairing is expected to shift 
to 3 P 2 pairing. We treat the suppression of the specific heat 
following Maxwell (1979). For core neutrino processes (direct 
Urea, modified Urea, and nucleon bremsstrahlung), we multi- 
ply the normal state neutrino emissivity by one Boltzmann 
factor exp {-A/kT), where A is the relevant gap, for each par- 
ticipating superfluid nucleon. 

Neutron ‘So pairing occurs in the inner crust only. Crust 
neutron superfluidity has little effect on the cooling of a young 
neutron star. It may become important if older neutron stars 
are reheated by the motion of vortices unpinned in pulsar 
glitches (Alpar et al. 1984), but we do not include this effect in 
our model. Thus, our calculations will underestimate the 
surface temperatures of older neutron stars. We take our treat- 
ment of l * S 0 neutron pairing from Ainsworth, Wambach, & 
Pines (1989). 

Proton *S 0 pairing is the subject of some controversy. 
Several authors (Chao, Clark, & Yang 1972; Takatsuka 1973; 
Niskanen & Sauls 1981; Amundsen & 0stgaard 1985a) have 
found that proton superfluidity occurs, but that the S 0 gap 
vanishes when the proton Fermi momentum reaches ~ 1.0— 12 
fin -3 . If this result is true, proton l S 0 pairing will have almost 
no effect on direct Urea cooling because any equation of state 
which allows the direct Urea process to operate will have 
proton Fermi momenta much larger than this. A possibly dif- 
ferent result is indicated by the work of Wambach, Ainsworth, 
& Pines (1991). These authors include a treatment of polariza- 
tion effects and find that the proton *S 0 gap is smaller than 
that found in previous calculations, but that the gap does not 
seem to decrease with increasing density. Because of this uncer- 
tainty we have performed a calculation using the gap found by 
Chao et aL (1972), which has almost no effect, and calculations 
with several density independent transition temperatures, 
which have a large effect because the superfluidity extends 
throughout the direct Urea core in the inner part of the star. 

Neutron *P 2 pairing has been calculated by Hoffberg et al. 
(1970), Takatsuka (1972), and Amundsen & 0stgaard (1985b). 
The transition temperatures found in these calculations differ 
by almost an order of magnitude. In addition, Takatsuka 
(1972) finds that the gap vanishes at a neutron Fermi momen- 
tum of 2 fin -1 while the other two calculations find that the 
gap persists to a much higher density. The chief culprit in the 
uncertainty seems to be the treatment of dispersion effects. In 
addition, none of the calculations has included a treatment of 
polarization effects, and preliminary calculations by Jackson et 
al. (1982) indicate that they strongly enhance the pairing. We 
present calculations using the neutron 3 P 2 gaps found in each 
of these calculations. 


3. RESULTS 

3.1. Observational Data 

There are several promising candidates for the detection of 
thermal X-rays from cooling neutron stars, but none of them is 
definitive. We plot the inferred surface temperature or the 
upper limit on the surface temperature for six of the best candi- 
dates on Figures 2 and 3 for comparison with our cooling 
curves. For Vela (Ogelman & Zimmermann 1989) and PSR 

1055 — 52 (Brinkmann & Ogelman 1987), a blackbody spec- 
trum gave the best fit to the EXOSAT data. The Einstein 


i i 



Observatory HRI failed to resolve the sources RCW 103 
(Tuohy et al. 1983) and PSR 0656+14 (Cordova et al. 1989), 
suggesting that the observed X-rays were emitted from the 
surface of the star rather than from a synchrotron nebula sur- 
rounding the star. However, ROSAT has not detected the 
compact source in RCW 103, giving thus an upper limit of 
\2 x 10 6 K on its temperature (Becker 1992). The surface tem- 
peratures for 3C 58 (Becker, Helfand, & Szymkowiak 1982) 
and the Crab (Hamden & Seward 1984) are upper limits based 
on Einstein data. 

All of the surface temperatures plotted in the figures are the 
temperature measured by an observer at infinity. 


3.2. Cooling without Superfluidity 
The dramatic effect of the direct Urea process on the thermal 
evolution of a neutron star is demonstrated in Figure 1, which 
shows the cooling curves for the neutron star models listed in 
Table 1. Nucleon superfluidity is neglected in all calculations. 
Our equation of state has a minimum mass of 135 Af© for the 
onset of the direct Urea process, so the models with Af £ 1.4 
Af© cool by direct Urea and the models with Af ^ 1.3 Af© do 
not The photon luminosity of a direct Urea star in the age 


TABLE 1 


Properties of Neutron Stars of Various Masses 



Rgtat 





Pc 

1.0 

11.910 

0.897 

9.76 

0 

0 

494 

1.1 

11.780 

1.007 

9.90 

0 

0 

348 

1 4 

11.630 

1.119 

9.99 

0 

0 

3.67 

1J 

11.460 

1.228 

10.03 

0 

0 

4.15 

14 

1 1-250 

1.339 

10.02 

0.038 

2.405 

4.76 

\S 

10.970 

1.448 

9.93 

0.212 

4.169 

5.58 

1.6 

10-570 

1.560 

9.72 

0.489 

5.345 

6.91 

1.7 

9.470 

1.674 

8.91 

1.003 

6438 

12.07 


Note.— Masses and radii of the stars, their cores, and, for Af > 1-35 Af©, 
their inner cores f P»0 where the direct Urea process is allowed. Masses are 
in solar and radii in km. The last column gives the central density in 

units of nuclear matter density p 0 ” 2.8 x 10 u g cm *. 
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range I0 2 -I0 5 yr is roughly I0 4 times smaller than a modified 
Urea star of the same age. Note that the direct Urea process 
has an enormous effect even if the fraction of the core in which 
it is allowed is very small (see the 1.4 M 0 case). The mass 
dependence of the surface temperature for the direct Urea stars 
of age 10 2 -T0 5 yr is due to the more massive stars having larger 
direct Urea cores. 

After a rapid initial cooling phase, which lasts roughly a 
year, the surface temperature stays on a plateau for several 
decades. The duration of the plateau phase is determined by 
the thermal time scale of the crust. The systematic effect of the 
more massive stars leaving the plateau the earliest is due to the 
fact that the more massive a neutron star is, the less massive a 
crust it has. The smaller crusts of the more massive stars have 
shorter thermal time scales. The level of the plateau is deter- 
mined by the properties of the crust, so the surface temperature 
during this phase of the cooling is independent of the cooling 
rate of the core. The dramatic difference between direct Urea 
cooling and modified Urea cooling appears once the crust has 
thermally relaxed and the star leaves the plateau. 

3.3. The Effect of Core Neutron Superfluidity 

Cooling curves for the 1.4 M 0 neutron star using various 
results for the 3 P 2 neutron pairing gap and allowing for the 
direct Urea process are given in the solid curves in Figure 2. All 
curves include crust neutron superfluidity with the gap 
calculated by Ainsworth et al. (1989) and proton superfluidity 
with the l S Q gap calculated by Takatsuka (1972), but neither of 
these gaps has any significant effect on the cooling. For com- 
parison, the dashed curve has no 3 P 2 neutron pairing and the 
two dotted curves correspond to a 1.3 M 0 star cooling without 
direct Urea, with and without 3 P 2 neutron pairing. 

The qualitative effect of core neutron superfluidity is to halt 
the rapid cooling of the interior once the interior has reached 
the superfluid transition temperature. This results in higher 



log t (years) 


Fig. Z — Effect of neutron superfluidity on the direct Urea cooling. Cooling 
curves for a 1.4 M 0 direct Urea star using various calculations of the neutron 
J Pj gap are shown in the solid curves. The labels are defined in the text. The 
dashed curve has no neutron y P 2 superfluidity. The dotted curves show the 
cooling of a 1.3 M 0 star without direct Urea, with and without neutron J P 2 
pairing for comparison. The sources of the plotted data are given in the text 



Log t (years) 

Fig. 3. — Effect of proton superfluidity on direct Urea cooling. Cooling 
curves for a 1.4 M 0 direct Urea star using various assumptions about the 
proton l S 0 gap are shown. The dash-dot curve uses the proton gap found by 
Chao et al. (1972). The dashed curve and the two solid curves have been 
computed using density-independent critical temperature for proton ‘Sq 
superfluidity as labeled. The dotted curves show the cooling of a 13 M 0 star 
without direct Urea, with and without proton pairing for comparison. 

surface temperatures once the crust has completed its thermal 
relaxation. Roughly speaking, the surface temperature stabil- 
izes at T t = TfxT e \ where 7; is the superfluid critical tem- 
perature, a is typically ^0.2, and TfT in J is the surface 
temperature-interior temperature relation found by Gud- 
mundsson, Pethick, & Epstein (1982, 1983) and Hemquist & 
Applegate (1984). 

This qualitative picture is borne out by the detailed calcu- 
lations presented in Figure 2. The cooling curves labeled 
HGRR, 0.3HGRR, and 0.1HGRR have the neutron 3 P 2 gap 
calculated by Hoffberg et al (1972), and that gap multiplied by 
0.3 and 0.1. The curve labeled AO uses the gap calculated by 
Amundsen & 0stgaard (1985b), and the curve labeled T72 uses 
the calculation of Takatsuka (1972). Core neutron superfluidity 
with the Takatsuka (1972) gap has almost no effect on the 
cooling because the gap vanishes in the direct Urea core of our 
star; there is no superfluidity where it is needed. For the rest of 
the models, the larger gaps give higher surface temperatures. 
Multiplying the Hoffberg et al. (1972) gap by 0.1 lowers the 
cooling curve by just about 0.5 in the log. The TfT ln J relation is 
7^ cc 7^ 55 , so the detailed result is exactly what the qualitative 
picture predicted. 

3.4. The Effect of Proton Superfluidity 

The effect of proton superfluidity on the cooling of the 1.4 
Af © star is illustrated in Figure 3. All curves have been calcu- 
lated using the l S 0 neutron gap of Ainsworth et al (1989) and 
the 3 P 2 neutron gap of Takatsuka (1972). Neither of these 
neutron gaps affects the cooling significantly. The dot-dash 
curve (CCY) has been calculated using the proton gap 
obtained by Chao et al (1972). This proton gap vanishes in the 
inner core in which the direct Urea process is allowed in this 
model, so proton superfluidity has almost no effect on the 
cooling, as can be seen by comparing with the triple dot-dash 
curve which has no proton superfluidity at all. The two dotted 
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curves are for comparison purposes. They correspond to a 1.3 
M 0 star, which cools without direct Urea, with and without 
proton superfluidity. The two solid curves give results of calcu- 
lations in which density-independent critical temperatures of 
1.2 and 2 x 10 9 K for proton superfluidity have been used, 
corresponding to the upper and lower values found by 
Wambach et al. (1991), and the dashed one has a higher critical 
temperature of 3 x 10 9 K. These curves show a large effect due 
to proton superfluidity because the superfluidity, by assump- 
tion, occurs throughout the direct Urea core in the center of 
the star. The one with T c = 3 x 10 9 K gives temperatures even 
larger than the models with only modified Urea processes and 
no superfluidity. 

4. CONCLUSIONS 

We have included the direct Urea process in calculations of 
the cooling of young neutron stars and found that stars in 
which the direct Urea process is allowed cool to invisibility as 
soon as the thermal relaxation of the crust has been completed 
unless nucleon superfluidity intervenes. This qualitative result 
is independent of the size of the inner core in which the direct 
Urea process occurs; direct Urea cooling is so efficient it domi- 
nates the cooling if it is allowed anywhere in the star. Super- 
fluidity of either neutrons or protons can shut off the direct 
Urea process and halt the drop of the surface temperature, but 
the entire direct Urea inner core must become superfluid. If 
superfluidity occurs and shuts off the direct Urea process, the 
surface temperature of the star stays essentially constant until 
the star is £ 10 5 yr old, at which time surface photon cooling 
takes over and the cooling speeds up. Direct Urea stars in the 
age range 10 2 -10 5 yr old are thermometers for nucleon super- 
fluidity because their surface temperatures are determined by 
the superfluid transition temperature in the direct Urea inner 
core and almost nothing else. 

Figures 2 and 3 show that the Vela pulsar and PSR 
0656 + 14 have surface temperature below the standard cooling 
curves, as already found by Nomoto & Tsuruta (1986, 1987) 
and Page & Baron (1990, 1991), and thus require some 
enhanced cooling controlled by nucleon superfluidity. Super- 
fluidity transition temperatures of the order of 10 9 K are then 


required down to the center of the star , if the reported tem- 
peratures of these two objects are indeed surface temperatures. 
Drawing a more quantitative conclusion will require the inclu- 
sion of nucleon isobars and hyperon processes and a more 
accurate treatment of the effect of *P 2 pairing on the neutrino 
emissivity, as well as better data. Older stars (e.g., PSR 
1055 — 52) may have significant heating from the motion of the 
superfluid vortex lines in their crusts. 

The nuclear symmetry energy is the most important factor in 
determining whether or not the direct Urea process is allowed 
in neutron stars. Our results have been obtained with an equa- 
tion of state whose symmetry energy was adjusted to give the 
onset of direct Urea cooling at 1.35 Af 0 . While this choice was 
made to illustrate the dramatic effect that direct Urea cooling 
can have, we do not believe that this choice seriously limits the 
significance of our results. Our neutron star models, sum- 
marized in Table 1, have quite standard radii and central den- 
sities. Our equation of state has a standard value for the 
compression modulus of symmetric nuclear matter, and our 
symmetry energy is consistent with that found by Lattimer et 
al. (1991) for relativistic nuclear equations of state. We con- 
clude, in agreement with Lattimer et al. (1991), that direct Urea 
cooling in neutron stars is consistent with all constraints on the 
nuclear equation of state, and its possible occurrence must be 
considered in neutron star cooling calculations. 

Neutron stars are probably born with a range of masses. If 
the range of iron core masses of their progenitors translates 
into a similar range of neutron star masses, then a mass range 
ofO.2-0.4 M 0 is expected. It is possible that the more massive 
stars cool by direct Urea and the less massive stars do not In 
this case neutron stars will divide into two populations, differ- 
ing substantially in their thermal histories. 

We thank T. L. Ainsworth, L. Bildstein, M. Prakash, & D. 
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